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In this work we examine the Fat-Link Irrelevant Clover (FLIC) overlap quark propagator and the 
gluon propagator on both dynamical and quenched lattices. The tadpole-improved Luscher-Weisz 
gauge action is used in both cases. The dynamical gauge fields use the FLIC fermion action for 
the sea quark contribution. We observe that the presence of sea quarks causes a suppression of 
the mass function, quark renormalization function and gluon dressing function in the infrared. The 
ultraviolet physics is unaffected. 
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Quark and gluon propagators are fundamental quanti- 
ties in QCD encoding the rich nonperturbative and per- 
turbative properties of QCD. In previous studies, the Fat- 
Link Irrelevant Clover (FLIC) overlap quark propagator 
and the duon propagator were examined on quenched 
lattices 3 . In this work we study the effects of two 
flavor dynamical FLIC fermions Q on these quantities 
for the first time. 

The use of smeared or "fat" links in lattice fermion 
actions has been of interest for some time 5] and FLIC 
fermions [6, 7] have shown a number of promising advan- 
tages over standard lattice actions. The FLIC fermion 
action is a Wilson-clover type fermion action in which 
the irrelevant operators of the Wilson and clover terms 
are constructed using fat links, while the relevant opera- 
tors use the untouched (thin) gauge links. By smearing 
only the irrelevant, higher dimensional terms in the ac- 
tion, and leaving the relevant dimension-four operators 
untouched, short distance quark and gluon interactions 
are retained. 

Fat links may be created by any number of smearing 
prescriptions including APE smearing [1, [^, [lo|. jlll , hy- 
percubic smearing [12], stout- link smearing [l3| or the 
unit-circle smearing .4;] used herein. 

The use of fat links in the FLIC action minimizes 
the effect of renormalization on the action improvement 
terms. Scaling studies indicate FLIC fermions provide 
a new form of nonperturbative 0{a) improvement [T^ 
where near-continuum results are obtained at finite lat- 
tice spacing. 

Access to the light quark mass regime is enabled by 
the improved chiral properties of the lattice fermion ac- 
tion [l5|. In particular, the histogram of the additive 
mass renormalization encountered in chiral-symmetry 
breaking Wilson-type fermion actions is seen to narrow 
upon introducing fat-links in the irrelevant operators. 
These benefits facilitate the generation of the dynami- 
cal fermion configurations examined herein. 



I. CONTINUUM PROPAGATORS 

We begin by reviewing the formulation of the quark 
and gluon propagator on the lattice. In the continuum, 
the tree- level {A^{x) = 0) quark propagator is identified 
with the (Euclidean space) fermionic Greens function. 



(1) 



where is the bare quark mass. In momentum space 
this equation is solved straightforwardly. 



Af(p) 



(2) 



Denote S^^\p) = Af(p) to be the tree-level propagator 
in momentum space. Then in the presence of gauge field 
interactions, define Sharcip) to be the fourier transform 
of the (interacting) fermionic Green's function. 



{]p + m°)A^^''%x, y) = S^{x - y). 



(3) 



We define the mass function M{p) and the bare renor- 
malization function Z{p) such that the bare quark prop- 
agator has the form 



5'barc(p) 



Z{P) 



ii) M(p) ■ 



(4) 



Then for the renormalization point C, the renormalised 
quark propagator is given by 



ii) + M (p) 



Z2iC,a)Shiii-c{p), 



(5) 



where Z,^{p) is the (C-dependent) renormalization func- 
tion, and ^2(^,0) is the wave function renormalization 
constant, which depends on ^ and the regulator param- 
eter a. The a-dependence of S'barc is implicit. Z2(C,a) is 
chosen such that 



1. 



(6) 
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As S^{p) is multiplicatively renormalisable, all of the (- 
dependence is contained within Zq{j)), that is, the mass 
function is (^-independent. 

The continuum tree-level gluon propagator in Landau 
gauge is associated with the following Greens function, 

{S^'d^ - d^^dnSatD^^^llix, y) = 5\x - y). (7) 

In momentum space the tree-level gluon propagator takes 
the form. 



(8) 



The non-perturbative gluon propagator is defined as the 
following two-point function. 



D';^,{x,y)^{Al{x)Al{y)). 



(9) 



The scalar propagator D{q^) is related to the full propa- 
gator (in momentum space) by 

<.('Z) = (^M."^)^"W)- (10) 

The renormalised scalar propagator D^{q^) is chosen such 
that at some momentum scale ^, 

q^D^{q^%.=^2 = 1. (11) 

We refer to q'^D^{q^) as the gluon dressing function. 

II. FLIC OVERLAP FERMIONS 

The overlap formalism [H, [13, HI, [ll] in the vector- 
like case leads to the following definition of the massless 
overlap-Dirac operator 20] , 



(12) 



Here, €{H) is the matrix sign function applied to the over- 
lap kernel H. The kernel can be any Hermitian version 
of the Dirac operator which represents a single fermion 
species of large negative mass. The standard choice is 
the Hermitian Wilson-Dirac operator (setting a = 1), 



i?w =75(y + ^A-mw), 



(13) 



where y is the central covariant finite difference operator, 
and A is the lattice Laplacian, or Wilson term. 

In this work, as in previous studies, we use the FLIC 
action [1,0] as the overlap kernel 1]. The FLIC fermion 
action is a variant of the clover action where the irrelevant 
operators are constructed using APE-smeared links [I,!, 
[Tol . [T]| . and mean field improvement [2l| is performed. 
The Hermitian FLIC operator is given by 

Hmc = l5{fmR + liAln - \<J ■ F^fi) - mw), (14) 



where the presence of fat (smeared) links and/or mean 
field improvement has been indicated by the super- and 
subscripts. We choose a^^ = f [7/ij Iv] and use a standard 
one- loop F^i,, 



ly^ — j-l 



{x) + U-^^-,{x)), 



(15) 



(16) 



where U^j,^y{x) is the elementary plaquettc in the 
direction. 

The APE smeared links Uj}{x) constructed from U^{x) 
by performing n smearing sweeps, where in each sweep 
we first perform an APE blocking step. 



V^/\x) = {l-a)^+^Y.^ 



(17) 



followed by a projection back into SU{3),ulf\x) = 

iv/l''\x)) ■ In this work, the projection is defined by 
first performing a projection into U{3), 



U'{V) = V[V^V]-i, 
and then projection into SU{3), 

1 



V{V) = 



{/det U'{V) 



U'iV). 



(18) 



(19) 



As it is only the product na that matters [24I, we fix 
a = 0.7 and only vary n. 

Mean field improvement is performed by making the 
replacements 



fl ' 



(20) 



where mq and Uq are the mean links for the standard and 
smeared gauge fields. We calculate the mean link via the 
fourth root of the average plaquette. 



1 i 
Mo = (-ReTr{7^^(a:)}^^,,<^. 



(21) 



III. LATTICE QUARK PROPAGATOR 

It is easily seen that the continuum massless quark 
propagator anti-commutes with 75, 



{75,5'{;arc(p)|„o=oJ' = ^■ 



(22) 



A straightforward consequence of the Ginsparg- Wilson 
relation is the inverse of the overlap Dirac operator sat- 
isfies 



{75,^0"'} = 275 



(23) 
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Noting that the mean hnk is a function of a and that 
Mo(a),Ug(a) ^ 1 as a ^ 0, we obtain 



hm Do 

a— >0 



2m„ 



(24) 



It is then natural to define the (external) massless bare 
overlap propagator on the lattice as 



!rnal) m 

Mm 



<S'baro(p)|mO=0 = 



1 



2mv 



as we then have that 

{75,'S'baro(p)|„o^o} = 0, 



(25) 



(26) 



as in the continuum case. The massive overlap Dirac 
operator is given by [l^l 



(27) 



with < 1 representing fermions of mass oc j-^- The 
massive (external) bare overlap propagator is defined 
as fH 



2mw(l — /i) 
and with the identification 



2mv 



satisfies 



(28) 



(29) 



(30) 



In order to construct M{p) and Z{p) on the lattice, we 
first define B{p),Cn{p) such that 



5barc(p)--«jC(p)+S(p). 



Then 



B{p) 



1 

UsUc 



■Tr[7^S'barc(p)], 



Tr[5^--(p)], 



(31) 

(32) 
(33) 



where the trace is over color and spinor indices only, and 
Us, ric specify the dimension of the spinor and color vector 
spaces. Now, define the functions B{p),Cf^{p) such that 
the inverse of the bare (lattice) quark propagator has the 
form 



s^Lip)^^(Pip) + B{p)- 

Then it is easily seen that 



C 

r = ^ 



B = 



B 



C2 + B^' 



(34) 



(35) 



where C"^ =C-C. 



The kinematical lattice momentum is defined such 
that at tree-level 



(5(°))-i(p) = z^ + mO, 



(36) 



that is qf,{p) = Cj°\p),m° = B^^^p). We note that the 
simple form of these relations is one of the advantages of 
overlap fermions, as the absence of additive mass renor- 
malization prevents the need for having to perform any 
tree-level correction [25j (necessary most other fermions 
actions [H, [13, [13)' outside of identifying the correct 
momentum variable q. Now, we define A{p) such that 



ShLip)^-^ 



B{p). 



(37) 



The mass function M (p) and renormalization function 
Z{p) may then be straightforwardly constructed, 



Zip) 



1 



A{p)' 



M(p) 



Bjp) 
A{p)- 



(38) 



IV. LATTICE GLUON PROPAGATOR 

We use a tadpole-improved plaquette plus rectangle 
gluon action. 



4aTgc = q/3E E I^^Tr[(l - C/^.(x))- 



\-,[2-R%^^\x)-R(^r\m- (39) 



20^2'" "'^ 



The lattice gauge field may be related to a continuum 
gauge field through a path ordered exponential. 



(40) 



We can recover the lattice gluon field through the follow- 
ing "mid-point" definition, 

l-TT{U^{x)~Ul{x)) + 0{a'). (41) 

As in the continuum, the lattice gluon propagator in 
coordinate space is then given by 



^;t(^,y) = (^^(^)4(y))- 



(42) 



A Fourier transform takes us to momentum space, and we 
calculate the scalar propagator D{q^) directly. Tree level 
improvement of the scalar gluon propagator is performed 
by identifying the appropriate kinematical momenta q 
such that 



i^(«)(g2)^ 1. 



(43) 



We denote the two different (gluonic and fermionic) kine- 
matical momentum by q, and use context to distinguish 
them. 
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V. SIMULATION DETAILS 



Calculations are performed on two quenched and two 
dynamical lattices. The dynamical lattices use a FLIC 
fermion action for the two degenerate flavors of sea quark. 
All the lattices have approximately the same physical vol- 
ume, and the details of all lattices are given in Table [H 
Two lattice volumes are used, 12^ x 24 and 16'^ x 32, with 
the lattice spacings for the quenched and dynamical lat- 
tices approximately matched. Landau gauge is chosen for 
the gauge fixing. An improved gauge fixing scheme [29| is 
used, and a Conjugate Gradient Fourier Acceleration |30| 
algorithm is chosen to perform the gauge fixing. 

Each lattice ensemble consists of 50 configurations. 
The lattice fermionic Green's function is obtained by in- 
verting the FLIC overlap Dirac operator on each config- 
uration using a multi-mass CG inverter [sH . We use pe- 
riodic boundary conditions in the spatial directions and 
anti-periodic in the time direction. The Fourier trans- 
form is taken to convert to momentum space, and the 
the bare quark propagator is obtained from the ensemble 
average. The quark propagator is calculated for 15 differ- 
ent masses. The details of the FLIC overlap parameters 
used are presented in Table [Til The matrix sign func- 
tion in the FLIC overlap is evaluated using the Zolotarev 
rational polynomial approximation [32], of degree (typi- 
cally 8) chosen to give an accuracy of 2.0 x 10^^ within 
the spectral range of the kernel, after projecting out low- 
lying modes. The tree level propagator is calculated by 
setting the links and the mean link to unity (the free the- 
ory). The kinematical lattice momentum q is obtained 
numerically from the tree level propagator, although it 
could equally well have been obtained from the analytic 
form for q derived in Ref. (25l |. 

The gluon propagator is obtained by performing a 
Fourier Transform of the lattice gluon field using the 
midpoint definition and then calculating the scalar prop- 
agator directly. 



Volume 


/3 




a (fm) 


rriTT 


uo 


Phys. Vol. (fm-*) 


12^ 


X 24 


4.60 


0.0000 


0.120 


CO 


0.8888 


1.44^ X 2.88 


16^ 


X 32 


4.80 


0.0000 


0.096 


oo 


0.8966 


1.54^ X 3.08 


12^ 


X 24 


4.00 


0.1318 


0.120 


806 


0.8338 


1.44^ X 2.88 


16^ 


X 32 


4.20 


0.1300 


0.096 


820 


0.8745 


1.54^ X 3.08 



TABLE I: Parameters for the different lattices. All use a tad- 
pole improved Luscher Weisz gluon action. Lattice spacings 
are set via a string tension analysis incorporating the lattice 
Coulomb term. Shown is the lattice volume, gauge coupling 
f}, lattice spacing, pion mass (for the dynaical quarks) in MeV, 
the mean link and the physical volume. 



VI. RESULTS 

We first examine the gluon propagator results. Fig- 
ure [T] shows the bare gluon dressing function q^D{q^) 
on each of the lattices. The squares represent momenta 
entirely in a spatial cartesian direction while triangles in- 
dicate points where the momentum is entirely along the 
temporal axis. The separation between these is due to 
the use of an asymmetric lattice and is a clear indication 
of finite volume effects. This is to be expected on such 
small physical volumes. The quenched and dynamical 
renormalised gluon dressing functions q^D^{q^) are com- 
pared directly in Figure [51 We see that the presence of 
dynamical quarks causes the suppression of the dressing 
function in the infrared. This effect is clear even with 
the relatively heavy sea quarks used in this study. The 
same effect was seen in work which compared the gluon 
propagator on quenched and dynamical lattices with a 
staggered fermion action[33l [33|. 

We now turn to the quark propagator results. The 
mass function M{p) and the renormalization function 
Zt(p) are calculated on each of the lattices. A cylinder 
cut [2^ is applied to all the data, to reduce the effects 
of rotational symmetry violation. The full results for the 
dynamical lattices are displayed in Fig. [3l for M{p) and 
Fig. [SI for (p) . The renormalization point for (p) 
is chosen to be approximately q = 6 GeV. The full re- 
sults for the quenched FLIC overlap quark propagator 
presented in previous work[l|. 

In order to compare the quenched and dynamical re- 
sults, for each fixed momenta value we perform a quar- 
tic fit of both M{p) and Z{p) as a function of the pion 
mass (excluding the lightest two masses for finite vol- 
ume reasons). The data and fits for the smallest 10 





a = 0.120 


a = 0.096 


Ml 


0.00400 


0.00305 


M2 


0.00800 


0.00610 


M3 


0.01200 


0.00915 


/i4 


0.01600 


0.01220 


fJ-s 


0.02000 


0.01525 


M6 


0.02400 


0.01830 


fJ-7 


0.02800 


0.02134 


MS 


0.03200 


0.02439 


^9 


0.04000 


0.03049 


Mio 


0.04800 


0.03659 


Mil 


0.06000 


0.04574 


M12 


0.08000 


0.06098 


M13 


0.10000 


0.07623 


M14 


0.12000 


0.09148 


M15 


0.14000 


0.10672 



TABLE II: The 15 lattice mass parameters jj, used for each a. 
Values are chosen so that the bare masses are approximately 
matched for each lattice. 



5 



momenta on each of the lattices are shown in Fig. [4] 
for M{p) and Fig. [5] for Z{p). Having obtained the 
fit fmictions M(p, rn^) and Z{p, mV) we then choose 
ml 0.0, 0.25, 0.5, f. 0,2.0, 3.0 GeV^ and compare the 
functions at these values of the pion mass. In this way 
we can compare the quenched and dynamical results with 
matched pion masses. Results are shown in Figure [71 to 
which we now turn our attention. 

First we examine the mass function. The results for 
a = 0.120 show little effect due to the presence of sea 
quarks, just a slight suppression in the infrared for the 
lighter masses. The heavier masses show essentially no 
difference between the quenched and dynamical results, 
execpt in the extreme ultraviolet, where the discretisa- 
tion errors give some difference. The discretisation errors 
are not present in the a = 0.096 results, but now we can 
clearly see infrared suppression of M{p) in the dynamical 
results for all masses. 

Turning to the renormalization function, both the lat- 
tices show infrared suppresssion of {q) in the presence 
of dynamical quarks, and agreement in the far ultraviolet. 
At a = 0.096 the quenched and dynamical results also 
agree in the intermediate momentum regime, whereas 
there is a slight discrepancy in this region at a = 0.120. 

VII. CONCLUSIONS 

The gluon propagator and FLIC overlap quark prop- 
agator are compared on quenched and dynamical FLIC 



lattices at a = 0.096 and a = 0.120. We observe the 
the presence of dynamical quarks causes suppression in 
the infrared in all quantities examined: the gluon dress- 
ing function q^D{q^), the mass function M{p) and the 
quark renormalization function Z{q). Some differences 
between the results for the two lattice spacings indicate 
we are not quite in the scaling region at the coarser 
spacing. An alternative source of systematic error in 
the the quark propagator results is the simplistic fitting 
functions used, which do not take into account the pos- 
sibility of non-analytic behaviour at light quark mass. 
However, our results are consistent with those obtained 
elsewhereHlllllli]. 
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FIG. 1: Gluon propagator data for the bare gluon dressing function q^D{q^). The quenched lattices are shown on the top and 
the dynamical lattices on the bottom. The left column shows the results for a — 0.120 and the right column for a = 0.096. 
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FIG. 2; Plots comparing the quenched (squares) and dynamical (triangles) renormalised gluon dressing function, q^D^{q^), for 
a = 0.120 (left) and a — 0.096 (right). The renormalization point ^ was chosen to be 6.0 GeV. 



8 



1.0 

0.9 

0.8 

0.7 

% 0.6 

^ 0-5 h 
ft 

^0.4 
0.3 
0.2 
0.1 
0.0 



" * i i 



i Si i 



• I 

3 

n n 
I i 



i i i 5 1 i_ 

J J J i i i 



* i Si i i I i i 



5 5 5 5 



^ Iff 



2 3 
p GeV 




1.0 

0.9 

0.8 

0.7 

% 0.6 
O 

^ 0.5 
ft 

^0.4 



• • • «_ 



°^ O CO Q_ 



A A A A 



!!!!••• ** 



2 3 
p GeV 



FIG. 3: Cylinder cut data for the dynamical FLIC Overlap mass function M{p) at finite quark mass for the for the lattices at 
a = 0.120 (left) and a = 0.096 (right). The plots are against the discrete lattice momentum p. 




FIG. 4: Cylinder cut data showing the dependence of the mass function M(p) for the quenched (left) and dynamical (right) 
FLIC Overlap on the bare mass mP at fixed momenta, for the lowest 10 momenta values. At small bare mass, the curves are 
ordered inversely to the momenta they represent, that is the smallest momenta is the topmost curve. Data is shown for the 
lattices at a = 0.120 (top) and a — 0.096 (bottom). The solid curves are the quartic fits to the data. 
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FIG. 6: Cylinder cut data showing the dependence of the renormalization function Z{p) for the quenched (left) and dynamical 
(right) FLIC Overlap on the pion mass at fixed momenta, for the lowest 10 momenta values. The curves are ordered 
according to the momenta they represent, that is the largest momenta is the topmost curve. Data is shown for for the lattices 
at o = 0.120 (top) and a = 0.096 (bottom). The solid curves are the quartic fits to the data. 
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FIG. 7: Cylinder cut data comparing the interpolated mass function M(p) (left) and renormalization function Z^lq) (right), 
for the dynamical and quenched lattices at a = 0.120 (top), a — 0.096 (bottom). The different curves from lowest to highest 
represent matched pion masses of = 0.0, 0.25, 0.5, 1.0, 2.0, 3.0 GeV^. The solid points indicate the chiral limit. The different 
Z( (g) have been offset vertically for clarity. The renormalization point is chosen to be C « 6 GeV. 



